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We demonstrate that the NMSSM can have small fine-tuning and modest light stop mass while 
still evading all experimental constraints. For small tan/3 (large tan/3), the relevant scenarios are 
such that there is always (often) a SM-like Higgs boson that decays to two lighter — possibly much 
lighter — pseudoscalar Higgses. 



in 
o 
o 



> : 

in . 
o ■ 

<n : 
o . 
in ■ 

o ; 

^ : 

a: 
i 



X 



In the CP-conserving Minimal Supersymmetric Model 
(MSSM), large soft-supersymmetry-breaking mass pa- 
rameters are required in order that the one-loop correc- 
tions to the tree-level prediction for the lightest Higgs 
boson (m/j < mz) increase mn sufficiently to avoid con- 
flict with lower bounds from LEP data. The large size of 
these soft- SUSY breaking masses compared to the weak 
scale, the natural scale where supersymmetry is expected, 
is termed the little-hierarchy problem. This hierarchy 
implies that a substantial amount of fine-tuning of the 
MSSM soft-SUSY breaking parameters is needed. The 
severity of these problems has led to a variety of alter- 
native approaches. For instance, little Higgs models Q 
can be less fine tuned. Or, one can argue that large 
fine-tuning is not so bad, as in "split-supcrsymmctry" 
0. In this letter, we show that the Next to Minimal 
Supersymmetric Model (NMSSM can avoid or at 
least ameliorate the fine-tuning and little hierarchy prob- 
lems. In addition, we find that parameter choices that 
are consistent with all LEP constraints and that yield 
small fine-tuning at small tan /? (large tan (3) are nearly 
always (often) such that there is a relatively light SM-like 
CP-even Higgs boson that decays into two light, perhaps 
very light, pseudoscalars. Such decays dramatically com- 
plicate the Tevatron and LHC searches for Higgs bosons. 

The NMSSM is very attractive in its own right. It 
provides a very elegant solution to the [i problem of 
the MSSM via the introduction of a singlet superfield 
S. For the simplest possible scale invariant form of the 
superpotential, the scalar component of S naturally ac- 
quires a vacuum expectation value of the order of the 
SUSY breaking scale, giving rise to a value of [i of order 
the electroweak scale. The NMSSM is the simplest su- 
persymmetric extension of the standard model in which 
the electroweak scale originates from the SUSY break- 
ing scale only. A possible cosmological domain wall 
problem £| can be avoided by introducing suitable non- 
renormalizable operators [j| that do notgenerate dan- 
gerously large singlet tadpole diagrams Hence, the 
phenomenology of the NMSSM deserves to be studied at 
least as fully and precisely as that of the MSSM. 

Radiative corrections to the Higgs masses have been 
computed 0, 0, 0, 0] and basic phenomenology of the 
model has been studied [Hj . The NMHDECAY program 
|l2j allows easy exploration of Higgs phenomenology in 



the NMSSM. In particular, it allows for the possibility 
of Higgs to Higgs pair decay modes (first emphasized in 
(l3j and studied later in Q ) and includes the associated 
modifications of LEP limits. Of greatest relevance are 
h — > aa decays, where ft is a SM-like CP-even Higgs 
boson and a is a (mostly singlet) CP-odd Higgs boson. 
The relevant limits come from the analysis llfil of the 
Zh — > Zaa — > Zbbbb channel and the analysis |l6( of the 
t~ channel. The weaker nature 



Zh — > Zaa — > Zt + t 



of the limits from LEP on such scenarios will play an 
important role in what follows. 

The extent to which there is a no-lose theorem for 
NMSSM Higgs discovery at the LHC has arisen as an 
important topic 0, 0, 0, 0, H3 ■ In particular, it has 
been found that the Higgs to Higgs pair decay modes can 
render inadequate the usual MSSM Higgs search modes 
that give rise to a no-lose theorem for MSSM Higgs dis- 
covery at the LHC. And, it is by no means proven that 
the Higgs to Higgs pair modes are directly observable at 
the LHC, although there is some hope [l8lll^. 

Earlier discussions of fine-tuning in the NMSSM have 
been given in 0, H3- 

We very briefly review the NMSSM. Its particle con- 
tent differs from the MSSM by the addition of one CP- 
even and one CP-odd state in the neutral Higgs sector 
(assuming CP conservation), and one additional neu- 
tralino. We will follow the conventions of [l^. Apart 
from the usual quark and lepton Yukawa couplings, the 
scale invariant superpotential is 



(1) 



A SH u H d + — S 3 
depending on two dimcnsionless couplings A, k beyond 
the MSSM. [Hatted (unhatted) capital letters denote su- 
pcrficlds (scalar superfield components).] The associated 
trilincar soft terms are 

XA x SH u H d + ^A K S 3 . (2) 

The final two input parameters are 

tan/3 = h u /hd , Mcff = As , (3) 
where h u = (H u ), h d = (H d ) and s = (S). These, 
along with mz, can be viewed as determining the three 
SUSY breaking masses squared for H u , H d and S (de- 
noted Trijj, m 2 H and rn/g) through the three minimiza- 
tion equations of the scalar potential. 

Thus, as compared to the three independent param- 
eters needed in the MSSM context (often chosen as /Lt, 
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tan/3 and M A ), the Higgs sector of the NMSSM is de- 
scribed by the six parameters 

A , k , A x , A K , tan/3 , ^ off . (4) 
We will choose sign conventions for the fields such that 
A and tan/3 are positive, while k, A\ 1 A k and // ff should 
be allowed to have either sign. In addition, values must 
be input for the gaugino masses and for the soft terms 
related to the (third generation) squarks and sleptons 
that contribute to the radiative corrections in the Higgs 
sector and to the Higgs decay widths. 

Sample discussions of the fine-tuning issues for the 
MSSM appear in [23|. We will define 



F = Max a F a = Max a 



d\ogm z 



dloga 



(5) 
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where the parameters a comprise /x, and the other 
GUT-scale soft-SUSY-brcaking parameters. (In some pa- 
pers, d ^" n z is employed.) In our approach, we choose 
m^-scalc values for all the squark soft masses squared, 
the gaugino masses, M 1j2) 3(toz), A t (mz) and A^imz) 
(with no requirement of universality at the GUT scale). 
We also choose m^-scale values for tan/3, n and rriA', 
these uniquely determine B^mz). The vevs h u and hd at 
scale mz are fixed by tan/3 and mz via m 2 z = ~g 2 (hu + hd) 
(where g 2 = g 2 + g' 2 ). Finally, m 2 H ^(rnz) and m 2 Hd {mz) 
are determined by the two potential minimization condi- 
tions. We then evolve all parameters to the MSSM GUT 
scale (including fi and B^). Next, we shift each of the 
GUT-scale parameters in turn, evolve back down to scale 
mz, and reminimizc the Higgs potential using the shifted 



FIG. 1: Left: the fine-tuning measure F in the MSSM is 
plotted vs. Im~ m~ , without regard to LEP constraints on 

V J l J 2 

uih- The + points have ran < 114 GeV and are excluded 
by LEP limits. The x points have mj, > 114 GeV and are 
experimentally allowed. Right: F is plotted vs. nrih for all 
scanned points. 

terms of /x 2 , tan/3, m 2 H ^ and m 2 H , with \x replaced by /x e ff- 
However, a second combination gives an expression for 
Heff in terms of m 2 z and other Higgs potential parameters. 
Eliminating /j e ff, we arrive at an equation of the form 



,+1Bm 7 -\-C = 0, with solution ; 



-B±y/B 2 - C, 



where B and C are given in terms of the soft susy break- 
ing parameters, A, k and tan (3. Only one of the solutions 
to the quadratic equation applies for any given set of pa- 
rameter choices. Small fine-tuning is typically achieved 



values of /Lt, £? M , m 2 Hi and mj Id . This gives new values when C«B 2 and derivatives of m 2 z with respect to a 

GUT scale parameter tend to cancel between the —B and 
+\/B 2 - C (-VB 2 - C) for B > (for B < 0). 

To explore fine-tuning, we proceed analogously to the 
manner described for the MSSM. At scale mz, we fix 
tan/3 and scan over values of A < 0.5 (A < 0.7 is re- 
quired for perturbativity up to the GUT scale), \k\ < 
0.3, sign( K ) = ± and 100 GeV < \fx eS \ < 1.5 TeV, 
sign(/i e ff) = ±- We also choose TO^-scale values for the 
soft-SUSY-breaking parameters A\, A K , A t = Ab, Mi, 



for h u and hd yielding new values for mz and tan /3. 

Results will be presented for tan (3(mz) = 10, 
M lt2t3 (m z ) = 100,200,300 GeV. We scan randomly 
over |j4((m^)| < 500 GeV and 3rd generation squark 
and slepton soft masses-squared above (200 GeV) 2 , as 
well as over > 100 GeV, sign(/i) = ± and over 

to A > 120 GeV (for which LEP, MSSM constraints re- 
quire m h > 114 GeV H3). On the left side of Fig. [TJ 
we plot F as a function of the mean stop mass ^J' 



m—^rn 



which enters into the computation (we use HDECAY |2 
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with m\ ole = 175 GeV) of the radiative correction to the 
SM-likc light Higgs mass mh- Points plotted as +'s (x's) 
have mh < 114 GeV (mh > 114 GeV) and are excluded 
(allowed) by LEP data. Very modest values of F (of 
order F ~ 5) are possible for mh < 114 GeV but the 
smallest F value found for mh > 114 GeV is of order 
F ~ 185 [j^. The very rapid increase of the smallest 
achievable F with mh is illustrated in the right plot of 
Fig. ^ This is the essence of the current fine-tuning 
problem for the CP-conserving MSSM. Also, to achieve 



i2, TOq, m\j, mjj, m\ , and m 2 E , all of which enter 
into the evolution equations. We process each such choice 
through NMHDECAY (using to^ ' = 175 GeV) to check 
that the scenario satisfies all theoretical and available ex- 
perimental constraints (including m~ > 100 GeV). For 
accepted cases, we then evolve to determine the GUT- 
scale values of all the above parameters. The fine-tuning 
derivative for each parameter is determined by shift- 
ing the GUT-scale value for that parameter by a small 
amount, evolving all parameters back down to mz, rede- 
termining the potential minimum (which gives new val- 



m h > 114 GeV, 



> 1 TeV is required, an indi- ucs h' u and h' d ) and finally computing a new value for 



cator of the little hierarchy problem. 

We now contrast this to the NMSSM situation. One 
combination of the three potential minimization equa- 
tions yields the usual MSSM-likc expression for m\ in 



using to^ 2 



9 2 (K 2 + h' d 2 ). 
Our results for tan/3 = 10 and Mi t 2,i(mz) = 
100,200,300 GeV and randomly chosen values for the 
soft-SUSY-brcaking parameters listed earlier are dis- 
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tan/S=10, M 123 (m z ) = 100,200,300 GeV 



250 



200 - 



150 




100 



500 750 1000 
(mrm r2 ) 1/2 (GeV) 



1250 1500 



FIG. 2: For the NMSSM, we plot the fine-tuning measure 

F vs. /m~ m~ for NMHDECAY-accepted scenarios with 
V ti i 2 ^ 

tan/3 = 10 and Mi, 2 , 3 (m z ) = 100,200,300 GeV. Points 
marked by '+' ('x') escape LEP exclusion primarily due to 
dominance of hi —> aiai decays (due to irih 1 > 114 GeV). 



tan/S=10, M 12i3 (m z ) = 100,200,300 GeV 
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FIG. 3: For the NMSSM, we plot the fine-tuning measure F 
vs. m/i-L for NMHDECAY-accepted scenarios with tan/3 = 10 
and M w (mz) = 100,200,300 GeV. Point labeling as in 
Fig.0 ' 



played in Fig. [2 We see that F as small as F ~ 5.5 can 
be achieved for m~ ~ 250-^400 GeV. In the figure, 
the + points have mi ll < 114 GeV and escape LEP exclu- 
sion by virtue of the dominance of hi ~ > a\a\ decays; as 
noted earlier, LEP is less sensitive to this channel as com- 
pared to the traditional hi — > bb decays. Points marked 
by x have > 114 GeV and will escape LEP exclu- 
sion regardless of the dominant decay mode. For most of 
these latter points hi — > bb decays are dominant, even if 
somewhat suppressed; hi — > aia\ decays dominate for a 
few. For both classes of points, the hi has fairly SM-likc 
couplings. We also note that all points with F < 20 have 
m hl < 114 GeV and BR(hi — > a±ai) > 0.70. Finally, in 
Fig. [HI we demonstrate the rapid increase of the minimum 
F with mj r The lowest F values are only achieved for 
Wfii Si 105 GeV. However, even for > 114 GeV, the 
lowest F value of F ~ 24 is far below that attainable for 



m h > 114 GeV in the MSSM. 

A small value for A K (mz) (typically of order a few 
GeV) appears to be essential to achieve small F. First, 
small A K allows small enough m ai that hi — > aiai 
decays are dominant; this makes it possible for the natu- 
rally less fine-tuned values of < 114 GeV to be LEP- 
allowed. Second, small F is frequently (nearly always) 
achieved for m/ H < 114 GeV (m^ > 114 GeV) via the 
cancellation mechanism noted earlier, where C <C B 2 , 
and this mechanism generally works mainly for small A K . 
Indeed, there are many phenomcnologically acceptable 
parameter choices with m\ lx > 114 GeV that have large 
A K , but these all also have very large F. 

For lower tan/3 values such as tan/3 = 3, extremely 
large Jmj m~ is required for ra^ > 114 GeV in 
the MSSM, leading to extremely large F. Results in 
the NMSSM for tan/3 = 3 arc plotted in Fig. g| for 
^1,2,3(7712) = 100,200,300 GeV and scanning as in the 
tan/3 = 10 case. We see that F ~ 15 is achievable for 
300 GeV. No points with m hl > 114 GeV 



m~ ~ 300 GeV. No points with m hl > 
were found. All the plotted points escape LEP limits 
because of the dominance of the hi — * aia\ decay. For 
very large tan/3 (e.g. tan/3 ~ 50), it is possible to ob- 
tain to/, > 114 GeV with relatively small ^/m~ m~ in 
the MSSM as well as in the NMSSM. We have not'yet 
studied fine-tuning at very large tan /3 in either model. 



tan|8=3, M 123 (m z )= 100,200,300 GeV 
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FIG. 4 

Fvs. y 



For the NMSSM, we plot the fine-tuning measure 
n~ for NMHDECAY-accepted scenarios with 

tan/3 = 3 and Mi, 2 , 3 (m z ) = 100,200,300 GeV. Point label- 
ing as in Fig. [5] 

In the NMSSM context, the smallest achievable value 
for F is mainly sensitive to M^(mz). For example, for 
M 3 (m z ) ~ 700 GeV and tan/3 = 10, the smallest F we 
find is of order F ~ 40. 

We note that in [2l| the mass of the SM-like Higgs 
h (where h = hi for the parameter choices they fo- 
cus on) is increased beyond the LEP limit by choosing 
modest tan/3 ~ 2 -j- 5 and A values close to the 0.7 up- 
per limit consistent with perturbativity up to the GUT 
scale. This maximizes the additional NMSSM tree-level 
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contribution to m\ proportional to A 2 , thereby allowing 
mi > 114 GeV for somewhat smaller , /to— to— than in 

V Cl 12 

the MSSM. This, in turn, reduces the fine-tuning and 
little hierarchy problems, but not nearly to the extent 
achieved by our parameter choices. In our plots, the SM- 
like h is always the h\. The points with very small F 
have low ^to~ to— , modest A and n, and escape LEP 
constraints not because vrih is large but because h — > aa 
decays are dominant. 

In conclusion, we reemphasize that the NMSSM pro- 
vides a rather simple escape from the large fine-tuning 
and (little) hierarchy problems characteristic of the CP- 
conserving MSSM. However, the relevant NMSSM mod- 
els imply a high probability for hi — > a\a\ decays to be 
dominant. We speculate that similar results will emerge 
in many supersymmetric models where the Higgs sector 
is more complicated than that of the MSSM. Higgs detec- 
tion in such a decay mode should be pursued with greatly 
increased vigor. Existing work 0, 0] which suggests a 
very marginal LHC signal for WW — > h\ — > a\a\ — > 
bbr + T~ when m ai > 2mi, should be cither refuted or im- 
proved upon. In addition, the a\a\ — > t + t~t + t~ chan- 
nel that dominates for 2m T < to Qi < 2to& (an entirely 
acceptable and rather frequently occurring mass range 
in our parameter scans and not excluded by T decays 
since the a% has a large singlet component) should re- 
ceive immediate attention. Hopefully, we will not have 
to wait for Higgs discovery at an e + e~ linear collider via 
the inclusive Zh — > £ + £~X reconstructed Mx approach 
(which allows Higgs discovery independent of the Higgs 
decay mode) or at a CLIC-based 77 collider |2(| in the 
77 — > h — > bbr + T~ or t + t~t + t~ modes. 
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